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ABSTRACT 

We investigated the thermal evolution of rotating strange stars with the de- 
confinement heating due to magnetic braking. We consider the stars consisting 
of either normal quark matter or color-flavor-locked phase. Combining decon- 
finement heating with magnetic field decay, we find that the thermal evolution 
curves are identical to pulsar data. 

Subject headings: Stars:strange star 
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Introduction 



Combining observational data available from the successful launch of Chandra and 
XMM/Newton X-ray space missions with the significant progress in theoretical studies of 
the physics of dense matter, these developments offer the hope for distinguishing various 
competing neutron star thermal evolution models. Through this way we would have 
an opportunity to determine various important properties of dense matter, such as the 
composition, superfluidity and the equation of state. These progresses should help us obtain 
deeper insight into the properties of dense matter. 

Quantum chromodynamics(QCD) predicted the existence of quark matter at 
high density. Moreover phenomenological and microscopic studies also confirmed 
that quark matter at a sufficiently high density, as in compact stars, undergoes 
a phase transition into a color superconductivity state, which are typical cases of 
the 2-flavor color superconductivity (2SC) and color-flavor locked (CFL) phases. 
Theoretical approaches concur that the superconductivity order parameter, which 
determines the gap in the quark spectrum, lies between 1 and lOOMev for baryon 
densities existing in the interiors of compact stars. It is generally believed that 
appearance of quark matter would be implied in structure and evolution of compact 



stars. Mass-radius relati on and changes in s t ructure have 



an amount of literatures 



Zheng fc Yu 20 06b 



Alcock 1986 



jeen deeply investigated in 



Glendenning et al. 1995a 



Glendenning et al. 1995b 



Yang fc Zheng 2002 



Pan fc Zheng20071 ). The evolutionary 



prop erties associated with spin o: 



work ( Bildst en fc Ushomirskv 2 000 



Pan et al. 20061 : IZheng fc Liu 2004 



;he star have been extensively discussed in past 



Madsen 2000 



Andersson 2000 



Zheng fc Kang 2005 



Zheng fc Pan 2006 



Zheng et al. 20071 ). Meanwhile 



someones had ever ex pected compact star cooling as a prob e to distinguish quark matter 



from hadronic matter tlLiu fcZheng 2005 



Zheng fc Yu 2006al ). The researches showed that 
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the cooler stars contain perhaps quark matter. However, the problem is unsolvable. As 
known, continuous deco nfinement processes occur during the rotatio n evolution of the 



star if quark core exists (I Zdunik et al. 2001 



Alford 2004! : 



Kang 20071 ). The corresponding 



deconfinement heating(DH) extremely influences the evolution temperature of the star. No 
matter whether a strange star which sustain a tiny nuclear crust is in normal phase or in color 
superconducing phase, DH should delay the co oling of the star. Consequently, the str ange 



star become warmer so that it wouldn't cooler fjYu &: Zheng 2006 



Zheng fc Zhou 2009 1 



It is well know that a compact star spins down due to magnetic dipole radiation. So the 
DH processes are closely related with the magnetic field. When we take constant field into 
account for the calculations of heating rate, the star would maintain too high temperature 
at old age evolutionary time, during which we have not found such hot pulsar yet (see figures 



in (lYu fc Zheng 2006 



Zheng fc Zhou 20061 )) except a few of millisecond pulsars. However, 



the observations support the fact that the field could decay such as Ohmic decay and the 
decay time scale is about 10 6 years. Obviously, the DH will decrease due to field decay at 
old ages. In this paper, we will assess how the thermal evolution of a strange star would be 
affected. 



2. Neutrino Emissivities and Specific Heat 



The most efficient cooling process in unpaired quark matter is the quark direct 
Urac(QDU) process d — > uev and ue — > di>, given by (jlwamoto 19821 ) 



8.8 x l0 26 a c uY^ /3 T!>( D ergon'* sec' 1 . 



where a c is the strong coupling constant, u = pbj Po,Pb is the baryon density and 

Po = 0.17 fm~ 3 is the nuclear saturation density,y e = p e /pb is the electron fraction, and T 9 is 

the temperature in units of 10 9 K. When the QDU process being switched off due to a small 
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electron fraction (Y e < Y ec = (3/7r) 1 / 2 mga c 3 ^ 2 /64.), the dominating contribution to the 
emissivities is the quark modified Urc a(QMU) dq — > uqev and quark bremsstrahlung(QB) 



Q1Q2 — ► Q1Q2W processes, estimated as ljlwamoto 19821 ) 



e (M) ~ 2.83 x 10 19 'a 2 c uT*( M ergcm^sec- 1 , 



e {QB) ^ 2.98 x 10 19 uT*( QB ergcm^sec- 1 . 



(2) 
(3) 



Because of the pairing in color superconducting phase, the emissivity of QDU process is 
suppressed by a factor of (d ~ exp(—A/T) and the em issivities of QMU and QB processes ar e 



suppressed by a factor (m ~ exp(—2A/T) for T < T ^Blaschke et al. 200 



)0; 



Shovkovv 2004J ). 



In order to compute the cooling cu rves of the stars, we need to give the specific heat of 



the electrons and quarks (jlwamoto 1982h : 



c e ~ 2.5 x 10 2 V /3 T 9 ergcm^R- 1 
c q ~ 0.6 x 10 20 y e 2/3 M 2/3 T 9 ergcm^R- 1 



(4) 



(5) 



But in c olor superconductivit y phase, the quark specific heat is changed 
exponentially (lBlaschke et al. 20001 ) 

2" 



c sg — 3.2c g ^ — I 



2.5 - 1.7 



3.6 



exp 



A 



(6) 



where T c is critical temperature related to A as A = 1.76T C . Compare with the total 



mass of the stars, the mass of the crust is very sma 



contribution to neutrino emissivity and specific heat (]Lattimer et al. 1994] ) 



\(M r < 1Q- 5 M^. So we neglect its 
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3. DH with Magnetic-field decay 



We consider a strange star embodied by a nuclear crust. The DH is determined by the 
mass change of the crust. The total heat released per unit time as a function of t is: 



tldec{t) — —Qn : — V, 



(7) 



nib dv 

where q n , the heat release per absorbed neutron, is expected to be in the range 
g n ~10 — 40MeV. Its specific value depending on the assumed SQM model, and m& is the 
mass of baryon. The mass of the crust M r can be approximated by a quadra tic function of 



rotation frequency v. As discussed in ( jZdunik et al. 2001 
the crust reads 



Yu fc Zheng 20061 ) . the mass of 

(8) 



M c = M c °(l + 0.24^ 2 + 0.16z/|) 
where u 3 = z//10 3 Hz and M c ° < 1O" 5 M is the mass of the crust in the static case. 



Assuming the spin-down is ind uced by the magne tic dipole radiation, the evolution of 



the rotation frequency v is given by ljYu &: Zheng 20061 ) 



2tt 2 



B(t) 2 R e v 3 sm 2 9, 



3/c 3 



(9) 



where / is the stellar moment of inertia, 9 is the inclination angle between magnetic and 
rotational axes. In our work, we combined the heating with magnetic field decay, so in 
Eq. ([nD the magnetic is denoted as a function of time. 

The field decay of the magnetic field is expected to be a very compli cated process. Fo r 
purposes of illustrate its main features, we assume a simple model equation (|Jose et al. 20071 ). 

dB B 

— = 10 

dt r D v 1 

with the initial condition B(0) = Bq. Eq. fllOIU s a crude approximation, we are dealing with 
a typical value of the magnetic field, neglecting the spatial variation of the field throughout 
the crust. Moreover, we adopt that the magnetic field dissipate on a time scale td 
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4. Cooling Curves 



Considering the energy run away and heating effect of the star, the cooling equation 
can be written as: 



(11) 



where Cy is the total specific heat, L u is the total neutrino luminosity and L 1 is the surface 
photon luminosity given by 



L 7 = 4tt j R 2 cxT s 4 , 



(12) 



where a is the Stefan-Boltzmann constant and T s is the surface temperature. The last term 
in Eq. (fT2l represents the DH due to the spin-down of the star. 

The surface temperature of the stars is related to internal temperature by a coefficient 
determined by the s cattering processes occurri ng in the crust. We apply an formula which 
is demonstrated by (IGudmundsson et al. 19831 ). It reads 



T s = 3.08 x lO^^Tg - 5495 , 



(13) 



where g Si i4 is the proper surface gravity of the star in units of 10 14 cm s 2 (IPotekhin et al. 200ll ) 



n principle, magneti c fields may change the expression of Eq ffTB"]) . However, 



Potekhin et al. 200 ll ) have present that the effect is negligible if the field strength 



is lower than 10 13 G. So Eq ffT3]) is a good approximation for our case. 

Considering the gravitational red-shift, and then the effective surface temperature 
detected by a distant observer is 



T°° = T, 



1-0.295 



M 



Ra 1 



lV2 



(14) 



here Rf, is the radium of the star in units of 10 6 cm 



We consider a canonical strange star of 1.4M at a constant dens ity in our wo rk, which 



is a very good approximation for strange stars of mass M < 1.4M0 (jAlcock 19861 ). We 
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choose u — 3, q n = 20 Mev, the initial temperature T = 10 9 K, initial period Pq = 0.78 ms, 
initial mass of the crust M c ° = 1O _5 M , the magnetic tilt angle 6 = 45° and the time scale 
r D = 2 x 10 6 yr. 

Using the model described in the preceding section, we plot the cooling curves 
of rotating strange star in normal phase with DH for various initial magnetic f ields 
(10 11 — 10 13 G) in Fig.l. The observational data are taken from (IPage et al. 2004 ). We take 



Y" e = 10 5 for Y e > Y ec which is a representative for the QDU process contributing to the 
cooling, whereas Y e = for Y e < Y ec when QMU and QB processes dominates. 

We also show the therm al evolution curves o f strange star in CFL phase in Fig.2. 
In contrast to previous work tjBlaschke et al. 20001 ) . DH increases surface temperature of 
strange stars effectively. With constant fields, the stars still maintain high temperature at 
old ages but we can't find such hot pulsars during the phase. In the case of field decay, the 
cooling curves are exponentially suppressed at tails. This modification is very important. 
Under consideration, the thermal evolutionary model is compatible with the pulsar data. 



5. Conclusion 

We have given out the cooling curves of rotating strange stars by considering the 
DH effect with magnetic field decay. The thermal evolution of rotating strange stars are 
different from previous works. The DH could increase the temperature of strange star to the 
inferred pulsar data while the magnetic field decay suppresses DH to retain rapid cooling at 
old ages(t > 10 6 yr). Obviously, there is no evidence for the existence of extra hot source 
at old age. 

We should noticed that the used model of magnetic field decay which is crude. We 
here apply a simple magnetic field decay model. The decay time scale is taken according 
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to statistical data of pulsar. Actually, the magnetic field decay model is an endless 
controversial problem. The rigorous discusses may be necessary but our conclusion won't 
be changed by future study. 



This work was supported by the NFSC under Grant No. 10603002. 



-10- 
REFERENCES 

Alcock, C, Farhi, E. k Olinto, A., 1986, ApJ, 310, 261. 

Alford, M. 2004, J. Phys. G, 30, 441. 

Andersson, N., 2000, ApJ. 534, L75. 

Bildsten, L., k Ushomirsky, G., 2000, ApJ. 529, L33. 

Blaschke, D., k Klahn, T., k Voskresensky, D.N. 2000, ApJ, 533, 406. 

Glendenning, N. K., Kettner, C, k Webber, F., 1995, Phys. Rev. Lett., 74, 3519. 

Glendenning, N. K., Kettner C, k Webber F., 1995, ApJ, 450, 253. 

Gudmundsson, E. H., Pethick, C. J., k Epstein, R. I. 1983, ApJ, 272, 286. 

Iwamoto, N., 1982, Ann. Phys., 141, 1. 

Jose A, P., Bennett, L., Juan A.M., k Ulrich, G., 2007 , Phys. Rev. Lett., 98, 071101. 
Kang, M., 2007, MNRAS, 375, 1503. 

Lattimer, J. M., Van Riper, K. A., Prakash, M. k Prakash, M. 1994, ApJ, 425, 802. 
Liu, X. W., k Zheng, X. P., 2005, Astropart.Phys., 24, 92. 
Madsen, J., 2000,Phys. Rev. Lett.. 85, 10. 

Page, D., Lattimer, J. M., Prakash, M. k Steiner, A. W., 2004, ApJS, 155. 
Pan, N. N., Zheng, X. P., k Yang S.H., 2006, MNRAS, 371,1359. 
Pan, N. N., k Zheng, X. P., 2007, Chinese J. Astron. Astrophys., 7,5L. 
Potekhin, A. Y., Yakovlev, D. G., k Prakash, M. 2001, A&A, 374, 213. 



Shovkovy, I. A., 2004, Lectures delivered at the IARD 2004 conference, Saas Fee, 

Switzerland, June 12-19, and at the Helmholtz International Summer School and 
Workshop on Hot points in Astrophysics and Cosmology, JINR, Dubna, Russia, 
Aug. 2-13. 

Yang, S. H., & Zheng, X. P., 2002, High Ener. Nucl. Phys., 26, 1172. 

Yu, Y. W., & Zheng, X. P., 2006, A&A, 450, 1071 

Zdunik, J. L., Haensel, P., & Gourgoulhon, E. 2001, A&A, 372, 535. 

Zheng, X. P., 2007, Chinese J. Astron. Astrophys., 6(Suppl.2), 129. 

Zheng, X.P, & Kang, M., 2005, Phys. Rev. C, 72, 025809. 

Zheng, X. P., & Liu, X. W., 2004, Phys. Rev. C, 70, 015803. 

Zheng, X. P., & Pan, N. N., 2006, New A, 12, 165. 

Zheng, X. P., & Yu, Y. W., 2006a, MNRAS, 369, 1071. 

Zheng, X.P, & Yu, Y.W., 2006b, A&A, 445, 627. 

Zheng, X. P. & Zhou, X., 2006, MNRAS, 371, 1659. 



This manuscript was prepared with the AAS L^TjtX macros v5.2. 



3.5 

1 2 3 4 5 
log(t[yr]) 



Fig. 1. — Cooling curves of rotating strange star in normal phase. The up group of lines 
are in the case which only mDURA process occured, the lower group of lines are in the case 
which DURA process switched on. The solid lines correspond to magnetic field decay effect. 
The dash-dotted lines represent the constant field case. 
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Fig. 2. — Cooling curves of rotating strange star in CFL phase. The up group of lines are 
in A = lMev, the lower group of lines are in A = O.lMev. The solid lines correspond to 
magnetic field decay effect. The dotted lines represent constant field case. 



